Commutability between Semiclassical Limit and Adiabatic Limit 
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We study the adiabatic limit and the semiclassical limit with a second-quantized two-mode model, 
which describes a many-boson interacting system. When its mean-field interaction is small, these 
two limits are commutable. However, when the interaction is strong and over a critical value, the 
two limits become incommutable. This change of commutability is associated with a topological 
change in the structure of the energy bands. These results reveal that nonlinear mean-field theories, 
such as Gross-Pitaevskii equations for Bose-Einstein condensates, can be invalid in the adiabatic 
limit. 
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The experimental creation of Bose-Einstein conden- 
sates (BECs) with dilute alkali atomic gases has gener- 
ated great excitement and literally created a new subfield 
in physicsUQ- One of the main reasons is that it has 
made it possible to test experimentally some fundamen- 
tal and important physics that could only be discussed 
theoretically before. For instance, Tonks-Girardeau gas 
and quantum phase transition between superfluid and 
Mott insulator has been studied by theorists since 1960s; 
they were observed experimentally only recently with 
BECs|i|i|. 

There are now even discussions on how to 
use BECs to study black holesQ and superstringsQ. 

In this Letter we discuss a fundamental concept in 
quantum mechanics, the commutability of the semiclassi- 
cal limit and the adiabatic limit, with a second-quantized 
two-mode model. We suggest a possible experimental 
testing of this concept with BECs. This concept is due 
to M.V. Berry 0. In brief, consider a quantum system 
whose Hamiltonian is time-dependent, 



ih 



dt 1 



H{R(at))\i>) 



(1) 



One can eliminate a from the above Schrodinger equa- 
tion with a scaled time t — at and an effective Plank 
constant h = ah. With this scaling argument, Hwang 
and Pechukas claimed that the semiclassical limit h — » 
and the adiabatic limit a — * are equivalent 

This point was refuted by Berry |7j, who pointed out 
that these two limits are not equivalent because the 
Hamiltonian H may depend implicitly on h. Moreover, 
he showed that these two limits are incommutable in a 
simple double-well model: the Landau-Zener (LZ) tun- 
neling rate in this model is zero if the adiabatic limit 
a — > is taken first; it becomes one when the semiclas- 
sical limit h — » is taken first 0- Since it is impossible 
to change h experimentally, this concept has remained a 
game of theorists. 

We revisit the commutability between the semiclassi- 
cal limit and the adiabatic limit with a second-quantized 
two-mode tunneling model. This model can be used to 



describe a BEC system where only two quantum states 
are important, such as in a double- well potential or with 
two internal quantum states 0, In this model, the 
semiclassical limit is N — > oo with N being the number of 
bosons. In this large N limit, the second-quantized model 
becomes a two-level mean-field model. We show that 
one can recover the second- quantized model by quantiz- 
ing this mean-field model with the Sommerfeld rule. As 
iV can be changed in experiments, the semiclassical limit 
becomes experimentally accessible. 

More interestingly, the commutability between the two 
limits, N — > oo and a — > 0, in this second-quantized 
model depends on its mean-field interaction strength c. If 
c is small, the two limits are commutable; when c is over a 
critical value, the two limits become incommutable. Such 
a dependence on c is found to be related to a topologi- 
cal change in the structure of the energy bands. These 
results indicate that nonlinear mean-field theories, such 
as Gross-Pitaevskii equations for BECs, can be invalid 
in the adiabatic limit when the mean-field interaction is 
strong. Finally, we discuss how this commutability can 
be tested in a BEC experiment. 

The second-quantized two-mode model is 



H 



7 



A 



(a f a - 6 f 6) + -{a)b + atf) - -{a^a - tfb) 2 , (2) 



where generators and annihilators a) , a and w , b are for 
two different quantum states. In the Hamiltonian H, 
7 is the energy offset between the two quantum states 
and changes with time as 7 = at. The parameter v 
measures the coupling between the two states while A > 
is the interacting strength between bosons. The minus 
sign before A indicates that the interaction is attractive. 
In this system the total number of bosons N is conserved. 

For this second-quantized model, its semiclassical limit 
is N — > 00. In such a limit, the system's dynamics is 
given by the following nonlinear two-level model, 



dt 



{[i-i<tf-|tr>. + h} 



(3) 



2 



where c = 7VA and \a\ 2 + \b\ 2 = 1. This model is of- 
ten called a mean-field model. Technically to obtain the 
mean-field model, one focuses on the Gross-Pitaevskii 
states^ \^ gp ) = "^r( a ^ + b&) N \v&c) . By computing 
the expectation value (H) = {^ gp \H\^ gp ), one obtains 
the mean-field Hamiltonian H m { — (H) /N (up to a triv- 
ial constant) in the limit of N — ► oo. The Hamiltonian 
H m { leads to the dynamics in Eq.©. For a rigorous ac- 
count of large N limit as a semiclassical limit in models 
such as Eq.J2J), we refer readers to Ref . [ll| . 

We emphasize that the semiclassical limit N — > oo is 
taken with the mean-field interaction strength c = NX 
kept constant. Physically, this is to ensure that the se- 
ries of systems with different N's have about the same 
physics. If A were kept constant instead of c, the last 
term in Eq. would become too dominating at the large 
TV limit, completely changing the physics of the system. 
When the model is used to describe a BEC in a 
double-well potential, the limit N — > oo at a constant 
c is equivalent to having a larger trap for more atoms 
in the BEC, or to tuning A smaller with the Feshbach 
resonance technique [l^. 

We are interested in how the second-quantized model 
Eq.© behaves in the two limits, N — > oo and a — > 0, 
in particular, whether the model's behavior depends on 
which limit is taken first. For this purpose, we follow 
Berry's methodology to focus on the tunneling behav- 
ior of the quantized model. 
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FIG. 1: Tunneling rate as a function of the mean-field inter- 
action strength c. The solid lines are obtained with the mean- 
field model © for a = 0.005, 0.001, 0.0001, 0.00001(from top 
to bottom); the circle-dashed lines are with the quantized 
model © for a = 0.005, 0.001, 0.0001(from top to bottom). 
N = 8 and v = 0.2 are used. 

In FigQ]the tunneling rates are plotted as a function 
of the mean-field interaction strength c. Two sets of 
tunneling rates are calculated: one with the quantized 



model |J2J for a fixed number of bosons; the other with 
the mean- field model ©. In computing the tunneling 
rate, we have assumed that the system is completely in 
state a at t — > —00; the tunneling rate is the probability 
of remaining in state a at t — * 00, the end of dynamical 
evolutions. At a fixed number of bosons, the dynamics 
of the quantized model can be found by expanding 
a quantum state in terms of Fock states \N a , Nb), where 
N a and Nb are numbers of particles in quantum states a 
and 6, respectively. In terms of \N a , Nb}, the Hamiltonian 
becomes a (N + 1) x (JV + 1) matrix. 

Upon careful examination of the data in Fig^ one no- 
tices that c = v is a critical value. When c < v, the 
tunneling rate goes to zero in the adiabatic limit a —> 
for both the mean-field model and the quantized model. 
However, when c > v, the tunneling rate from the mean- 
field model is always non-zero while the tunneling rate 
can be zero for the quantized model. Since the mean- 
field model is the semiclassical limit of the quantized 
model, the mean-field result can be regarded as the re- 
sult from the quantized model with the limit N — > 00 
having been taken. Therefore, the results in Fig^show 
that the tunneling behavior in the quantized model (0) 
depends strongly on the order of the limits taken while 
this dependence itself relies on the value of the mean-field 
interaction strength c. 




FIG. 2: Energy levels from the second-quantized model H 
(N — 20) and the mean-field model H m f. The solid lines are 
quantized energy levels; the open circles are mean-field energy 
levels. Note that for comparison with the mean-field theory, 
the quantized energy levels from H have been divided by TV. 

To understand the above results, we first examine the 
energy levels of the second-quantized model J2J as func- 
tions of 7, the slowly changing system parameter. These 
energy levels can be found by directly diagonalizing the 
Hamiltonian H and they are plotted in Fig|2 There is 
a drastic change in the structure of energy levels as the 
mean-field interaction c changes: a net of anti-crossings 
appearers in the lower part of the quantized energy levels 
when c > v. As known before [l^. when c > v there is a 
loop structure emerging in the energy band of the mean- 
field model J2J. When the mean-field energy levels (cir- 
cles) are also plotted in Fig|3 we find that the quantized 
energy levels are bounded by the mean-field energies. In 



3 



particular, the mean-field energy levels envelop the net 
of anti-crossings in the quantized energy levels. Such a 
correspondence was first noticed in Rcf. 9] . 

The structure change in the energy bands is associated 
with a change in the phase space of the mean-field model 
© as shown in Fig|3J In plotting this figure, we notice 
that the mean- field model, in fact, has only two indepen- 
dent variables and its Hamiltonian can be reduced to 



Hxal = 7P + o vi - 4p 2 cos ( 



cp 



(4) 



where p = (\a\ 2 — \b\ 2 )/2 and q = 6b — 6 a with 9 a j, being 
the phases of a and b. It is clear from FigOU when c < v, 
there is one minimum and one maximum; when c > v, 
we see two local minima, one maximum and one saddle 
point. Since these extremum points correspond to the 
eigenstates in the mean-field energy bands in Fig Etio], 
the structure change in the phase space is apparently 
connected with the structure change in the energy levels. 
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FIG. 3: Energy contours of the mean-field model _ff m f. Left: 
c = 0.1, v = 0.2, 7 = 0.0; right:c = 0.4, v = 0.2, 7 = 0.0. 

This connection can be further explored by re- 
quantizing the mean-field model H m f with the Sommcr- 
fcld theory, which says that the quantum motions are the 
periodic motions in the classical phase space that satisfy 

-!- / pdq = nh/N , n = 0, 1,2, ■•■ (5) 

2tt J 

The division by N comes from the fact that the mean- 
field Hamiltonian is an average for one particle, H m { = 
(H)/N. One can view K e ff = H/N as the effective Plank 
constant for H m {. In our calculations, the natural unit 
h = 1 is used. For convenience, we shall call the energy 
levels obtained with Eq.JSJ the Sommerfeld energy levels. 
They are shown and compared to the quantized energy 
levels of H in FigEI 

When c < v, the mean- field Hamiltonian has exactly 
one maximum (q = 0) and one minimum (q = tt). The 
Sommerfeld quantization around the maximum produces 
energy levels lower than the maximum energy while the 



FIG. 4: Comparison between the energy levels of the second- 
quantized model (dashed line) with N — 40 and the Sommer- 
feld energy levels (open circles), (a) c = 0.1, v — 0.2; (b) 
c = 0.4, v = 0.2. For clarity, we have plotted only a portion 
of the energy levels. 



quantization around the minimum generates energy lev- 
els higher than the minimum. This explains why the 
mean-field energy levels bound the quantized energy lev- 
els in Fig|2| We also see that the energy gap arises from 
the different quantization number in Eq.Q, from which 
we estimate that the energy gap between the lowest two 
energy levels at 7 = is A w vWl — c/v, independent of 
N. This agrees well with the numerical results in Fig[SJ 

When c > v, the phase space of H m f becomes very 
different: there are two local minima with an additional 
saddle point. In this case, the Sommerfeld quantization 
around the two local minima gives arise to two sets of 
Sommerfeld energy levels. In the lower part of Fig. 
for clarity, we have plotted only one set. If two sets were 
plotted, they would form a net of crossings, matching 
very well with the anti-crossing net from H. In doing 
the Sommerfeld quantization, we have ignored the tun- 
neling through the energy barrier between the two local 
minima. Once the tunneling is considered, degeneracies 
are lifted and the crossings become anti-crossings. This 
shows the energy gaps inside the triangular net have a 
different origin from the energy gaps outside the net or 
in the case of c < v. The energy gaps produced at these 
crossings can be estimated with the WKB method. Since 
the effective Planck constant for H m { is H/N, we expect 
that the gaps decrease exponentially with N. This is 
exactly what the numerical results in FigEI indicates. 

It is now not difficult to understand the tunneling be- 
havior that we have seen in Fig^ Let us recall the 
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FIG. 5: Energy gap between the lowest two eigen-energies in 
the second-quantized model at 7 = 0. The squares are for 
c = 0.1 and the dots for c = 0.4, with v = 0.2 for both. T he 
solid line is an approximation result A = uyl — c 2 /v 2 for 
c < v. 



Landau-Zener tunneling in a two-level model^jjj]. As 7 
changes with time as 7 = at, the LZ tunneling rate is 
Hz = cxp ( — , where A is the energy gap between 

the two levels. For a multi-level system like our second- 
quantized model H, the above equation should still be a 
very good approximation for the tunneling rate between 
two consecutive energy levels. We use the tunneling be- 
tween the two lowest energy levels as an example. As 
already analyzed, the energy gap changes with N as fol- 
lows, 



A = 



C < V . 



Nk2 exp(— T]N) c > v . 



(6) 



The parameter k\ « y/l — c/v. The other two parame- 
ters K2 and ry can be computed with the WKB method as 
in Ref.0 or with a more sophisticated method[l5|. This 
leads to the following tunneling rate 



r ~ r lz 



/ TTKiv \ 

/ ttN 2 k% \ 
eXp l~2^vJ C>V - 

For the case of c < v, it is clear that we have 



lim lim r = lim lim r = , 

N^ooa^Q a— >0JV— >oo 



(7) 



(8) 



which shows that the two limits a — > and N — > 00 are 
commutable. This explains why when c < v, both sets of 
the tunneling rates in FigQ become zero as a — > 0. 

For the other case c > v, the tunneling rate takes dif- 
ferent values at two different limits: 



lim lim r = , 
lim lim r > . 

a^0 N^co 



(9) 



This reveals that the two limits are no longer com- 
mutable. In the first limit, the adiabatic limit a — > 



is taken at a fixed number of bosons, for which the en- 
ergy gap is finite and one can always be slow enough 
not to causing tunneling. In the second limit, since the 
energy gap is already closed at N — > 00, tunneling oc- 
curs no matter how slow 7 changes. This explains why 
the tunneling rate from the mean-field model is always 
non-zero for c > v. This incommutability of these two 
limits also implies that the mean-field theories, such as 
Gross-Pitaevskii equation for BECs, can be invalid for 
the adiabatic limit. One example is the Bloch states for 
a BEC in an optical lattice is studied. In such a system, 
the Bloch wavenumber k can be regarded as an adiabatic 
parameter. If the Gross-Pitaevskii equation were always 
valid in the adiabatic limit, it would mean that stable 
Bloch states should exist for all possible k. However, as 
shown in Ref.^(|, a significant portion of Bloch states 
are unstable. 
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FIG. 6: Ratio of the bosons in the right well at the end of 
tunneling process. The computation is done with the second- 
quantized model H with a sweeping rate a — 0.0001. 



With cold atomic gases, we believe that it is now possi- 
ble to test experimentally the commutability between the 
semiclassical limit and the adiabatic limit. For instance, 
one can load a BEC into a double-well potential gener- 
ated by two laser beams 0. The energy offset 7 can 
be created by using different intensity for the two laser 
beams. To keep the mean-field interaction parameter c 
constant for different boson numbers, one can either use 
the Feshbach resonance technique |l^ to adjust the inter- 
action between atoms or change the size of the trap. In 
experiments, it is hard to measure the tunneling rate r 
between the two lowest energy levels as we just discussed. 
However, one can easily measure the number of atoms in 
either of the two potential wells. Once the experiment is 
set up, the most striking observation will be as shown in 
FigHj] For c < v, if the system initially has all its atoms 
in the left well (quantum state a) , then all the atoms will 
be in the right well for a fixed but very small sweep rate. 
It does not depend on N. For c > v, not all the atoms 
will fall into the right well: larger ./V less atoms in the 
right well. This difference illustrate our theoretical re- 
sults on the commutability of the two limits, a — > and 
N -> 00. 
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